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The relevance of anisotropic interactions in colloidal systems has recently emerged in the context 
of rational design of novel soft materials. Theoretical studies have predicted the possibility of a gas- 
liquid phase separation confined at low densities and the formation of empty liquids and equilibrium 
gels in low- valence systems. Here we provide experimental evidence of this scenario in Laponite, a 
complex colloidal clay with discotic shape and anisotropic interactions. We also report simulations 
of a patchy model for Laponite platelets, able to reproduce the observed experimental phase diagram 
and structural properties, confirming the crucial role of the reduced valence. 



Patchy colloids [TJ [2] are the novel building blocks of 
a bottom-up approach toward a rational design of self- 
assembled materials. The novelty of these particles re- 
sides in their intrinsic anisotropy, implemented either in 
shape or in the inter-particle interaction. Recent develop- 
ments in the synthesis of colloids of different shapes, pat- 
terns and functionalities [1 suggest that spontaneously 
self-assembled bulk materials with pre-defined properties 
can become feasible in the near future [3H3 • The possibil- 
ity of tuning the interaction anisotropy will make possible 
to recreate the molecular world at the nano- and micro- 
scale (a case with tremendous technological applications) 
as well as to generate novel unconventional phases, both 
ordered and disordered. Recent theoretical studies [8 , 
supported by numerical simulations [9], suggest, for ex- 
ample, that the phase-diagram can be significantly al- 
tered by limiting the particle coordination number (va- 
lence). New concepts such as empty liquids [8] (i.e. liquid 
states with vanishing density) and equilibrium gels [9- 
[llj (i.e. arrested networks of bonded particles, which do 
not require an underlying phase separation to form [12]) 
have emerged. So far, no experimental evidence of these 
predictions, specific of limited- valence systems, has been 
provided. 

In this work we describe the first experimental realiza- 
tion of the empty liquid scenario in a colloidal suspension. 
We study Laponite, an industrial synthetic clay made of 
nanometer-sized discotic platelets with positive charges 
on the plates and negative ones on the rims. The par- 
ticle shape and the inhomogeneous charge distribution 
generate a highly anisotropic interaction potential be- 
tween the discs. Similarly to other colloidal clays [131 |T4], 
Laponite has technological applications in cleansers, sur- 
face coatings, ceramic glazes, personal care and cosmetic 
products, including shampoos and sunscreens p^. It is 
characterized by a very rich, but controversial p!5UT9] , 
phase diagram including disordered (gels and glasses) 
and ordered (nematic) phases, on varying colloidal vol- 



ume fraction, at fixed ionic strength. The most debated 
region is the low concentration one (weight concentra- 
tion C^ < 2.0%), where the system ages very slowly up 
to a final non-ergodic state [181 US- This low concen- 
tration region is further investigated in this work, ex- 
tending the observation time to time-scales significantly 
longer than the previous ones. We discover that, despite 
samples appear to be arrested on the second time scale 
(as probed by dynamic light scattering [I8l|l9]), a signifi- 
cant evolution takes place on the year timescale. Samples 
indeed undergo an extremely slow, but clear phase sep- 
aration process into clay-rich and clay-poor phases that 
are the colloidal analog of gas-liquid phase separation. 
Spectacularly the phase separation terminates at a finite 
but very low clay concentration, above which the sam- 
ples remain in a homogeneous arrested state. These fea- 
tures of the phase diagram are strikingly similar to those 
predicted in simple models of patchy particles [8 , sug- 
gesting that Laponite forms an (arrested) empty liquid 
at very low concentrations. Differently from gels gener- 
ated by depletion interactions [T2j [20] or from molecular 
glass- formers [ST, where arrest occurs after the phase 
separation process has generated high-density fluctua- 
tion regions, here phase separation takes place in a sam- 
ple which is already a gel. Through a very slow rear- 
rangement, the system can thus explore the bottom of 
its free-energy landscape, finally reaching a stable min- 
imum which corresponds to an equilibrium gel, i.e. a 
disordered solid. 

Fig. [1JA.-C shows photographs of the temporal evolu- 
tion of a low concentration Laponite sample {Cuj=OA%). 
The initially fluid suspension (waiting time ^^^=0) 
(Fig. flK) progressively ages, forming a gel (the sam- 
ple does not flow if turned upside down as evident from 
Fig. [l^). The gelation time depends on clay concentra- 
tion and it is of the order of few thousand hours for low 
concentration samples [18]. Waiting significantly longer 
time (several years), the sample undergoes a thermody- 
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FIG. 1: Experimental behavior of diluted Laponite suspensions: (A-C) Photographs of a Cw — 0.4% sample (A) in the initial 
fluid phase (tw — 0), (B) in the gel state (tw ^ 4000 h), (C) in the phase-separated state (tw ~ 30000 h). (D) Photographs of 
samples in the concentration range 0.1 < Cw < 1.2% at very long waiting times (about 30000 h). All samples with Cw < 1.0% 
show a clear evidence of two coexisting phases, separated by an interface whose height (dashed horizontal lines) increases 
progressively with Cw (E) Estimated concentration of the denser gel phase in the separated samples shown in panel D. (F) 
Evolution of the measured S^ (Q) with waiting time for a Cw = 0.8% sample (located inside the phase separation region). (G) 
Evolution of S^ (Q) with waiting time for a Cw = 1.2% sample (located outside the phase separation region). The curves in 



panels F and G are measured at increasing waiting times. From bottom to top: tu 
8700, 11000 h. 



500, 900, 1600, 2700, 3400, 4700, 6000, 



namic phase separation, creating a sharp interface be- 
tween an upper transparent fluid and a lower opaque 
gel (Fig. [T]C). Phase separation is observed for all sam- 
ples with C^ < 1.0%. Fig. [ijl) shows a photograph of 
different concentration samples about three years after 
their preparation. The height of the colloid-rich part (in- 
dicated by the dashed lines in Fig. [Ijl)) increases pro- 
gressively with Cw, filling up the whole sample when 
Cw ^ 1.0%. This value thus marks the threshold of 
the phase separation region. Indeed, for concentrations 
larger than C^ = 1.0%, samples do not show any phase 
separation but maintain their arrested and homogeneous 
character at all times (see highest C^ sample in Fig. [Ijl)). 

An estimate of the concentration of the colloid-rich 
phase is provided by the ratio between the nominal con- 
centration and the volume occupied by the dense phase, 
i.e. neglecting the gas phase concentration. The esti- 
mated gel concentration C^^^\ reported in Fig. [lE, is 
~ 1.0% for all the phase-separated samples, as expected 
for phases which have reached their thermodynamic equi- 
librium concentrations. Furthermore, no macroscopic 
changes are observed, in the entire concentration range, 
in the following four years (seven years, around 60000 
hours, in total), strongly suggesting that all samples have 
reached their long-time equilibrium structure. The ab- 



sence of further restructuring, in samples above the co- 
existing liquid density, is consistent with the proposed 
equilibrium gel concept [8 . 

To characterize the evolution of the sample structure, 
we have monitored the static structure factor S^{Q) 
at different waiting times, from the initial fluid phase 
up to the gel state (arrested according to light scatter- 
ing measurements) and during the initial stages of the 
phase separation process. For this purpose we have per- 
formed Small Angle X-ray Scattering (SAXS) measure- 
ments [22 for more than one year on the same sample 
in the wavevector Q range 0.023 < Q < 0.4 nm~^ . We 
have focused on two different concentrations, respectively 
inside {Cw = 0.8%) and outside {Cw = 1.2% ) the phase 
separation region. The measured S^ {Q) are reported 
in Fig. [TJ^ and G for different waiting times tw- On in- 
creasing tw, S^{Q) of both samples show a shift of the 
main peak to higher Q values and an increase of the in- 
tensity at small wave vectors. However with the proceed- 
ing of the aging dynamics a drastically different behavior 
for the two samples is observed: while at Cw = 0.8% 
the small Q intensity continuously increases (Fig.[I^), at 
Cw = 1.2% the intensity saturates to a constant value 
at late times (Fig. [lb). These results suggest that the 
sample outside the coexistence region reaches its equilib- 
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FIG. 2: Behavior of the patchy particle model for Laponite discs. (A) Cartoon of a Laponite platelet and its schematization as 
a rigid disc composed by 19 sites (red spheres) with 5 attractive patches (blue spheres), three located on the rim and one at the 
center of each face. (B) Cartoon representing a T-bonded configuration for two interacting Laponite platelets and its realization 
in simulations. (C) Numerical phase diagram: binodal (blue curve) and percolation locus (red curve) in the p* — T* plane, 
where p* is the number density scaled by the close-packing density and T* is the thermal energy scaled by the strength of the 
bond (see Methods). (D-H) 3D snapshots of MC simulations at different state points. Different colors correspond to different 
clusters, and the red-color is reserved for the percolating cluster: D equilibrium fluid phase at T* = 0.10 and p* :^ 0.08; (E,F) 
equilibrium configuration at percolation for p* :^ 0.08 and 0.16; (G,H) final gel configurations at T* = 0.06 inside (p* :^ 0.08) 
and outside (p* :^ 0.16) the phase separation region. In these cases, platelets are connected into a single cluster (gel), which 
is clearly inhomogeneous (homogeneous) inside (outside) the binodal region. (I,L) Evolution of the S^'^^{Q) after a quench at 
T* = 0.06 for p* ~ 0.08, (inside the phase separation region) and p* ^ 0.16 (outside the phase separation region). Waiting 
times are: 10^ 1.2 x 10^ 5.7 x 10^ 1.6 x 10^ 3.6 x 10^ 6.1 x 10^ 10^ 2.2 x 10^ 4.9 x 10^ 1.1 x 10^ in MC steps. 



rium structure on the year-time scale, while the sample 
inside the unstable region, despite its apparent gel state, 
slowly evolves toward complete phase separation. Thus, 
in Laponite, gelation precedes, but does not preempt, 
phase separation. 

To favor the interpretation of the previous results we 
introduce a primitive model of patchy discs, which aims 
to mimic the strong rim-face charge attraction (signifi- 
cantly larger than the thermal energy) and the tendency 
to form open structures of Laponite clay. We aim at 
building a connection between the equilibrium phase di- 
agram of the model and the out-of-equilibrium dynam- 
ics and phase separation observed experimentally. Each 
Laponite platelet is schematized as a hard rigid disc com- 
posed by 19 sites on a hexagonal regular mesh, inspired 
by the work of ^^. To implement the rim- face link- 
ing (T-bonds [24) between different platelets, each disc 
is decorated with three sites on the rim and one at the 
center of each face (five sites in total). Only face-rim 
bonds can form, and they are modeled with a short-range 
square-well attraction, ensuring that each site can be in- 



volved at most in one T-bond. A representation of the 
model is provided in Fig.|2|^ and B. This primitive model 
highlights the anisotropic nature of the platelet-platelet 
interaction J23H25] but neglects the repulsive electrostatic 
barriers which control the timescale of the aggregation 
kinetics in Laponite. 

We perform Monte Carlo (MC) and Gibbs ensemble 
MC (GEMC) simulations to evaluate the gas-liquid co- 
existence region in the reduced density /?*- reduced tem- 
perature T* plane. Fig. [2p shows the binodal line, i.e. 
the locus of points separating homogeneous and phase- 
separated state, and the percolation line, defined as the 
line separating a finite cluster fluid phase (Fig. [2|l)) from 
configurations characterized by the presence of a span- 
ning infinite (transient) cluster. Fig. [2b and F show 
snapshots of the simulated system at the percolation line. 
As expected, the percolation line is located above the 
critical point J26\ and consistently with previously stud- 
ied patchy-spheres models [8 , the gas-liquid coexistence 
region is confined in a narrow window of T* and p* . In- 
deed, the coexisting liquid density, scaled by the closed 




FIG. 3: Phase diagram of diluted Laponite suspensions resulting from the combined experimental and numerical results in the 
concentration Cw and waiting time tw plane. Symbols correspond to experimental tw values requested to observe non ergodic 
behavior according to DLS |18|; boundaries inside colored regions are guides to the eye. For long waiting times, three different 
regions are identified, whose representative macroscopic behavior and a pictorial microscopic view are reported in (A-C) . (A) 
Phase-separated sample with colloid-poor (upper part) and colloid-rich (lower part) regions for Cw < 1.0%. The colloid rich 
phase is a realization of an empty liquid, arrested in a gel state. (B) Equilibrium gel for 1.0 < Cw < 2.0%, characterized by a 
spanning network of T-bonded discs. (C) Wigner glass, expected for 2.0 < Cw < 3.0% [30], where disconnected platelets are 
stabilized in a glass structure by the electrostatic repulsion, progressively hampering the formation of T-bonds. 



packed value, occurs at p* ~ 0.114, in significantly dilute 
conditions. Hence, it exists a wide region of densities, 
above the coexisting liquid density, where the system can 
be cooled down to very low T without encountering a 
phase separation, giving rise to an empty liquid state [8]. 
Such state is composed by an extensively-bonded per- 
colating network that, at low T, restructures itself on 
a timescale which exceeds the observation time, gener- 
ating an equilibrium gel state. We also study the out- 
of-equilibrium dynamics of the system. To mimic the 
experimental protocol, we first equilibrate the system at 
high T (corresponding to sample preparation) and then 
instantaneously quench it (corresponding to t^j = 0) to 
a sufficiently low T, so that the bond-energy is large as 
compared to the thermal energy. We propose to interpret 
the behavior of Laponite as the low-T limit of our model, 
connecting the increasing waiting time to a progressive 
temperature decrease in the numerical study [27l [28] . 
Snapshots of the final configurations of the system, af- 
ter a quench inside and outside the phase coexistence 
region, are shown in Fig. [2p and H. Independently from 
the density of the quench, the final configuration is always 
characterized by a single spanning cluster incorporating 
all particles. The structure of such cluster is highly in- 
homogeneous for quenches inside the coexistence region 
(Fig. [2b) and homogeneous for quenches in the empty 



liquid region (Fig. [2IH). Since at these low T the bond 
lifetime becomes much longer than the simulation time, 
the bonded network is persistent, i.e. the system forms 
a gel. 

Fig. [2J and L show the static structure factors calcu- 
lated from MC configurations S^'^^{Q) at several times 
(in MC-steps, equivalent to t^) following the quench, for 
two densities, respectively inside and outside the phase 
separation region. On increasing waiting time the scat- 
tered intensity increases in the region of the contact peak 
(T-bonds, Q ~ 0.5nm~^), revealing the bond formation 
process. The other notable feature is the increase of the 
scattering at small wavevectors. As in the experimen- 
tal data, two different scenarios occur at long times after 
preparation, respectively for samples inside (Fig. [2][) and 
outside (Fig. [2L) the unstable region. While inside the 
phase separation region S{Q) at small Q increases indef- 
initely, outside this region the growth stops after a finite 
waiting time, showing no further evolution. 

The zero-th order model introduced here for describ- 
ing Laponite at low densities condenses the electrostatic 
interactions between opposite charges into short-ranged 
attractive sites and neglects the overall repulsive electro- 
static interactions, following the spirit of primitive mod- 
els for charged molecules, an idea proposed long time ago 
when a patchy model for water was presented [29 . De- 



spite these simplifications, which reveal themselves in the 
different shape of S{Q) at short times — controlled by 
the screened electrostatic interactions ^i — and in the 
absolute values of S{Q) at small Q, the qualitative fea- 
tures shown by the model (Fig. 21 to L) coincide with 
the one measured experimentally (Fig. IF to G) both in- 
side and outside the unstable region. Most importantly, 
indefinite increase (saturation) in the growth of S{Q) at 
small Q is seen only inside (outside) the region where 
phase separation is observed, both in experiments and in 
simulations. 



The present results suggest a novel phase diagram 
of Laponite suspensions, summarized in Fig. [3) includ- 
ing the crossover taking place for C^ ~ 2.0% toward a 
Wigner glass, as recently suggested ^30j. The intricate 
free-energy landscape in the cross-over region between 
gel and glass has been investigated in Ref. [19]. At low 
concentrations, for C^ ^ 1.0%, the system evolves via a 
sequence (Fig. lA-C) of clustering (hours-days), gelation 
(months) and phase separation (years), from a sol to a 
homogeneous gel to a phase separated sample in which 
only the dense phase is arrested. This progression with 
tuj is strongly reminiscent of a constant-density path in 
the equilibrium phase diagram in which temperature is 
progressively decreased and the system evolves from a sol 
(FigJTK, Fig. |2t)), to a percolating structure (Fig. [lb, 
Fig.^E), finally encountering the phase-separation region 
(Fig7|l]C, Fig. |2]3). Differently from the case of isotropic 
short-range attractive colloids [12 where a homogeneous 
fluid is driven by a spinodal decomposition into an ar- 
rested network, here the system first forms a gel and then 
the gel extremely slowly increases its local density to ful- 
fill the search for a global free energy minimum (phase 
separated) state. These features are exactly the ones pre- 
dicted to take place in patchy colloidal systems when the 
average valence is small [27 , strongly supporting the view 
that the observed phase-separation is a genuine effect of 
the directional interactions. The gel phase observed in 
Laponite above 1.0% wt can thus be interpreted as an 
arrested empty liquid state, generated by the reduced va- 
lence, spontaneously arising from the combination of the 
platelet shape and the patchy distribution of opposite 
charges on the disc surface. The observed phase sepa- 
ration on the year time-scale calls attention on the fact 
that the long term stability of soft materials is controlled 
by the underlying phase diagram. Knowledge of thermo- 
dynamic properties is thus crucial in designing material 
with desired properties. Our case study shows that a 
careful choice of the density (within the empty liquid re- 
gion) may provide materials which are extremely stable 
in the long term (gels that do not phase separate nor age 
in the present case), since they are formed continuously 
from the liquid state, but finally reaching — through a 
very slow dynamics — their equilibrium configuration. 



Materials and Methods 
Laponite sample preparation 

Laponite RD suspensions were prepared in a glove box 
under N2 flux and were always kept in safe atmosphere 
to avoid samples degradation [32]. The powder, manu- 
factured by Rockwood Ltd, was firstly dried in an oven 
at T=400 C for 4 hours and it was then dispersed in pure 
deionized water {Cg — 10~^ M), stirred vigorously for 30 
minutes and filtered soon after through 0.45 jam pore size 
Millipore filters. The same identical protocol has been 
strictly followed for the preparation of each sample, fun- 
damental condition to obtain reliable and reproducible 
results [15]. The starting aging time {tw=^) is defined as 
the time when the suspension is filtered. 

Small Angle X-ray Scattering Small Angle X-ray 
Scattering (SAXS) measurements were performed at the 
High Brilliance beam line (ID2) at the European Syn- 
chrotron Radiation Facility (ESRF) in Grenoble, France, 
using a 10 m pinhole SAXS instrument. The incident x- 
ray energy was fixed at 12.6 keV. The form factor F{Q) 
was measured using a flow-through capillary cell. SAXS 
data were normalized and the scattering background of 
water was subtracted. The measured structure factor has 
been obtained as S^ {Q) = I{Q)/F{Q). 

Simulations 

Each platelet is modeled as a hexagon composed of 19 
hard-spheres of diameter a [23 , as schematically shown 
in Fig. |2|A.. A comparison with Laponite fixes <j = 5 nm. 
Each platelet is decorated with five sites, three located 
symmetrically on the rim and two on the two opposite 
faces of the central hard-sphere. 

Site-site interactions (acting only between rim and face 
sites) are modeled as square well interactions, with range 
0.1197<j and depth uq = 1. The number of sites controls 
the effective valence of the model. Since only rim-face 
bonds can be formed, the lowest energy state is charac- 
terized by an average number of bonds per particle equal 
to four. The choice of the valence controls the location 
of the gas-liquid unstable region, but does not affect the 
topology of the phase diagram [8 . Reduced temperature 
T* is scaled by the thermal energy, T* = /c^T/iio, where 
ks is the Boltzmann constant. Reduced density p* is de- 
fined as the number density p — N/L^ ^ where N is the 
number of particle and L the side of the cubic box, scaled 
by the closed packed density, corresponding to a hexago- 
nal close packing of discs (which is space-filling and equal 
to vT2)/19cr-3). Gibbs Ensemble MonteCarlo (GEMC) 
simulations are carried out for a system of 250 platelets 
which partition themselves into two boxes whose total 
volume is 66603<j^, corresponding to an average number 
density p* ^ 0.05. At the lowest studied T this cor- 
responds to roughly 235 particles in the liquid box and 
15 particles in the gas box (of side 32 a). On average, 
the code attempts one volume change every five particle- 
swap moves and 500 displacement moves. Each displace- 
ment move is composed of a simultaneous random trans- 



lation of the particle center (uniformly distributed be- 
tween zbO.OScr) and a rotation (with an angle uniformly 
distributed between ±0.1 radians) around a random axis. 

Standard MonteCarlo simulations (MC) are performed 
for a system of N =1000 platelets in the NVT ensemble. 
A MC step is defined as N attempted moves (defined as 
in the GCMC). Each state point is at first equilibrated 
at T* = 0.10, and then quenched down to T* = 0.06, a 
temperature well below the critical one, where the system 
cannot reach equilibrium within the duration of the run. 
The waiting time is defined as the time of the quench. 
To reduce numerical noise at each waiting time, the ob- 
servables of interest, such as the static structure factor 
S{Q) = 1/N{\pq\'^) with pq = e^^^, are averaged over 10 
independent runs. 

The use of a square-well potential to model the inter- 
actions makes it possible to unambiguously define two 
platelets as bonded when the pair-wise interaction en- 
ergy is —uq. Clusters are identified as groups of bonded 
platelets. To test for percolation, the simulation box is 
duplicated in all directions, and the ability of the largest 
cluster to span the replicated system is controlled. If the 
cluster in the simulation box does not connect with its 
copy in the duplicated system, then the configuration is 
assumed to be nonpercolating. The boundary between a 
percolating and a nonpercolating state point is then de- 
fined as the probability of observing infinite clusters in 
50% of the configurations. 
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